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:rj ■ Abstract 

5t \ The atomic parity-violation (APV) parameter Qt^ for a nucleus with n neutrons and z pro- 
tons has been included in the list of pseudo-observables accessible with the codes TOPAZO and 
ZFITTER. In this way one can add the APV results in the LEP EWWG 'global' electroweak 
fits, checking the corresponding effect when added to the existing precision measurements. 
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1 Introduction 



Recently we have been asked to include atomic-parity violation (hereafter APV) parameters 
in the list of pseudo-observables (hereafter PO) that are accessible with the FORTRAN codes 
TOPAZO im and ZFITTER [|], so to include the APV results in the LEP EWWG electroweak fits. 

The reason for this operation is that there are now precise experiments measuring APV in 
cesium 0, at the 0.4% level, thallium lead p and bismuth 0. Moreover, according to 0, 
the uncertainties associated with the atomic wave-functions have been reduced to another 0.4% 
for cesium. For additional uncertainties associated with the value of the tensor polarizability 
we refer to [Q. Note however that there is an intrinsic difference between the PO at the Z- 
resonance, e.g. F^, o"^, J^^^ etc, and the APV parameters where the typical scale is dictated 
by the limit of zero momentum transfer in the APV Hamiltonian. This fact alone is the origin 
of a comparatively larger theoretical uncertainty which is due to our basic ignorance of QCD 
corrections in this regime. 

The investigation of APV has been the subject of a number of studies made in the 80's 
by Marciano and Sirlin and |T^. For TOPAZO, which is based on the generalized minimal 



subtraction scheme |TT[, it has been relatively simple to include all recently computed higher- 



order effects in the old MS calculation. For ZFITTER instead, the authors have been able to 
produce a novel evaluation of the APV parameters in the on mass-shell (OMS) scheme. The 
current value for the weak-charge is 

Q^(Cs) = -72.06 ± 0.28 ± 0.34 (theo.) (1) 

For a recent evaluation of we refer, again, to where the program GAP? ||12| has been 
used. 



2 Upgrading the MS calculation 



The electron-quark parity-violating Hamiltonian at zero momentum-transfer will be conven- 
tionally parametrized as follows: 



where the ellipsis represents heavy-quark terms and we have factorized out the Fermi constant 
Gp. In heavy atoms the dominant part of parity-violation is proportional to the so-called 
weak-charge Qw- 

QUZ, A) = 2[{Z + A) Ciu + {2A-Z) Cid] , (3) 

We have taken the calculation by Marciano and Sirlin which is performed in the modified 
minimal subtraction scheme {MS) and have extended it to include all higher-order effects 
presently known. To summarize: two-loop leading contribution for the p-parameter , exact 



0{aag) corrections 
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aa"^ j corrections to p [|T§] , next-to-leading two-loop heavy top 
corrections [^. At the same time, an attempt has been made to evaluate the theoretical 
uncertainty at the level of electroweak and of QCD corrections. 

TOPAZO now returns, among all PO, the two quantities, Ciu and Cid of Eq.(0). They are 
defined as follows: 
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where sin^^(M^) is the MS weak-mixing angle at the scale = M^. We adopt a specific 
implementation of the re-summation procedure where the pair and sin^ 6'(M^) is the 
solution of a system of coupled non-linear equations that include all available higher-order 
effects, as described in Sect. 6.11 and 8 of ref. |T^. Moreover, 
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where = 1 — 4 s^, = sin^^(M^) and where we have two different treatments of the Z-7 
boxes — perturbative @ 
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and non-perturbative [|T2 
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where M is a mass scale representing the onset of the asymptotic behavior, i.e. the regime 
where becomes small. We observe a plateau of stability in K for M centered around 
0.5 GeV and this is the numerical value used. Furthermore we used the following form for the 
p, K, parameters : 
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In the previous equation we have h = M^/M'^. The strange quark mass is effectively chosen 
to be rus = 250 MeV so that (with effective rriu = md = 75 MeV) we recover the dispersive 
analysis for the Z-^ transition where n^-y(O) is rewritten in terms of a dispersion relation 
with the kernel connected to a{e~^e~ — > hadrons). Inside Eq.(|D 5ew(qcd) are the LO+NLO 
electroweak (^O (^a^ -|- QCD^ correction to p. The evaluation of p and includes the best 
available LO+NLO terms ^ . 

TOPAZO default is the perturbative formulation of the factorized result of Eq.@. There is 
the option of using some additive formulation where 
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where A^^ ^ are obtained from Eq.(^) by expanding and by neglecting terms of O (a^). 
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3 Atomic parity-violation in OMS scheme 



The old result of has been completely re-derived in the OMS scheme. Here, the technical 
problem is represented by the extraction of the limit of zero momentum transfer from the 
expressions that have been derived for the process ee ti Here the process under 

consideration is the t-channel scattering ee — ^ uu and what we need is naturally contained in 



the results of Ref. since they were derived retaining all masses and, therefore, the limit of 



zero momentum transfer, « (all) m? is possible. 

It is rather easy to take the limit for vertices and self-energy functions since they 

depend only on this variable. For boxes the procedure is more complex due to their complicated 
dependence on s and t invariants. Fortunately enough, ZZ and Z7 boxes form a gauge invariant 
sub-set of the whole result and for WW boxes one has to replace, in the corresponding limit, 
only the ^ = 1 part of the result, which is well defined and simple. This fact triggered the 
strategy for a calculation where we take all contributions but boxes from the — ^ limit of 
the ee ti form factors and were we have re-computed, from scratch, box diagrams at = 0. 
Note that this calculation was done with the aid of the computer system described in Ref. . 

For our calculation we compare the APV Hamiltonian of Eq.(^ with its ee ti analog, 
Eq.(I.lO) of 
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Here (0) stands for = and we write only one argument since box contributions are excluded. 
Moreover, 
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From Eqs.(|^) and (^) we immediately derive a relation between the APV parameters Cij and 
C2/ and the ee — > ti form factors at zero momentum transfer: 
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After a lengthy but straightforward calculations, we are able to reproduce the following generic 
expressions: 
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where we introduced a notation rij = mj/rrij and a fictitious term with non-zero neutrino 
charge in order to have a completely general representation and where 
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is a contribution, originating from the WW box, which is different for u and d channels (direct- 
crossed) . 

The other APV parameters can be obtained with the aid of some simple substitutions: 
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The terms of the second and third rows of Eq. ([1^) are identical to corresponding terms of the 
MS result. In the sequential order they are due to: QED vertex in Z exchange, W abelian 
vertex in 7 exchange (with neutrino charge), Z abelian vertex in 7 exchange; WW, Z7 and 
ZZ boxes. 

The only difference with respect to the MS result is present in the first term. The factor p 
is almost the same: 
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Ap'-(O) , (17) 



where we use instead the full expression for p (0), 

, Sf^^ (0) - J:^"" (0) , , 

Ap'''(0) = """^ ^ , (18) 

w 

contrary to the approximation made above where only the (leading) quadratic term in rrit is 
retained. The difference, being proportional to light fermion masses is numerically rather small. 

However, the main difference with the MS calculations is confined in the APV parameter 
Kpv for which, in the OMS-scheme, we derived: 
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where L^i^M"^ ) = In — ^. The gauge invariant Veltman Ap parameter is 
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and contains both the bosonic and the fermionic components. We explicitly give the bosonic 
part, Ap^°^ (for definition of finite part Bq of Bq functions see |]lj 
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To establish a link with the MS calculation we introduce the usual notion of leading and 
reminder terms: 
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where the leading term contains only Ap and the reminder contains all the rest: 
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Numerically, {AKpy)y^^^ and {AKpy)^.^^ are nearly equal and one might think that the usual 
leading-reminder splitting, the standard factorization of contributions with different scales and 
re-summation (see [0), 
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is not too well justified for the APV parameter Kpy. Note, however, that the factorized form 
of Eq.(p4D is fully consistent with the MS result Eq.(|) if we identify 



sin^ e^{M^ 



w 



(26) 



As done before, for n^J'^(O) we use effective quark masses which are consistent with a dispersive 



treatment of 11^^^ at zero scale. 

Finally, we apply mixed QCD (^O (^a^ + ct^)) and LO+NLO electroweak two-loop correc- 
tions for Veltman Ap parameter. For ppy we stick with one-loop (non-re-summed) result 
Eq.(|l^), since there the notion of leading-reminder splitting fails completely (numerically it 
looks like +3.5 — 3.0 = 0.5). For the same reason, we apply to Ap'^^''(0) only the mixed QCD 
but not the electroweak two- loop corrections, as already done in the first Eq.(^. The latter, as 
well as the other electroweak NLO corrections for remainder terms, although not implemented 
are successively used to evaluate the theoretical uncertainty in the electroweak sector of the 
OMS scheme. 



4 Theoretical uncertainty in APV 

In order to discuss the present level of theoretical uncertainty in atomic parity- violation we start 
with the MS results for Qiv(Cs) that are shown in Tab. |1|, corresponding to M^ = 91.1875 GeV, 
Mg = 150 GeV and ag(M^) = 0.119. ZFITTER numbers corresponding to 'Add/Pert' setup 
are added to the third row of the Table. The mj-dependence of Qyy{Cs) is shown in Tab. |l] 
where we register a 0.22(0.17) per-mill increase for between 170 GeV and 180 GeV and for 
Pert(Non Pert). As for the M^ dependence we have computed a decrease of about 0.7 per-mill 
for M„ between 150 GeV and 300 GeV. 
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Fact /Non Pert 


-73.1994 


-73.1932 


-73.1869 


Add/Pert 
ZFITTER 


-72.9732 
-72.9762 


-72.9658 
-72.9698 


-72.9582 
-72.9637 


Add/Non Pert 


-73.2026 


-73.1969 


-73.1912 



Table 1: Predictions for QwiG^i) from TOPAZO and ZFITTER for = 91.1875 GeV, = 
150 GeV and a^iMl) = 0.119. 



The associated theoretical uncertainty is approximately 3.2 per-mill and it is largely domi- 
nated by QCD effects. Let us consider the main sources of uncertainty. In the original calcu- 
lation of Marciano and Sirlin we have a dependence of the result on light quark masses. This 
appearance can be seen in Eq.(|^) and in the perturbative treatment of boxes, Eq.(|]). 

From 1983 the accuracy associated to the weak-charge Qw has been considerably reduced 
and we cannot include it in the list of high-precision PO if the result contains logarithmic 
enhancements due to light quark masses. 

In their second paper Marciano and Sirlin |10[ have suggested how to go beyond the partonic- 
language. One should distinguish quark masses in the Z-'j transition and in Z-7 boxes. Light 
quark masses, muds, are then fixed to parameterize the dispersive result for the Z-j transition 
and are not varied anymore in evaluating the theoretical uncertainty. 

Furthermore we have Z-7 box diagrams where quark masses show up as the consequence of 
the zero momentum transfer limit. Here, according to the suggestion of [T^ we split the boxes 
into a low-frequency part, approximated with the Born contribution for a physical nucleon 
(the term in Eq.(^), and an high-frequency part (the K-teim in Eq.(^) that includes 

O (ctg) corrections where light quark masses disappear. The mass scale M separating low- from 
high-frequency parts is, of course, arbitrary and only subjected to the requirement that a^lQ'^) 
starts to become small for | |> and that M > Aqcd- However, with the most complete 
evaluation of (up to three loops) we have found a plateau of stability for the result, i.e. for 
M between 0.5 and 0.6(0.8) K goes from 9.2016 to 9.1737(8.7818) and Q,„- has a variation of 
0.02(0.3) per-mill. Therefore we fix 0.5 < M < 0.6. 

Instead of varying light quark masses between undefined limits we prefer to estimate the 
theoretical uncertainty by comparing the perturbative result with light quark masses fixed to 
reproduce the dispersive approach to the Z-7 transition with a non-perturbative anzatz based 
on a low-frequency high frequency splitting at a mass scale of about 0.5 GeV. Note that when 
comparing B^p^M) with the perturbative factor In M^/M^ -|- | we find that the perturbative 
approach overestimates the effect of about 5.6%(1.9%) at M = 0.5(0.8) GeV. 

Furthermore, the differences in the factorized Eq.(^ versus additive Eq.(^ formulation of 
the coefficients Ciu^u is approximately 0.05 per-mill signaling that, from TOPAZO's treatment 
alone, pure electroweak higher orders are relatively under control. Another way of testing 
the electroweak theoretical uncertainty is, as usual, to compare two different renormalization 
schemes with the same input parameter set. When we compare ZFITTER in the preferred 
setup with TOPAZO additive/perturbative we obtain a relative difference of 0.04(0.05, 0.08) per- 
mill at rrit = 170(175, 180) GeV. However, an internal evaluation of electroweak theoretical 
uncertainties within ZFITTER (realized by evaluating the effect of the electroweak two-loop 
corrections which are not included in the preferred setup) shows a value of about ±0.25 per- 
mill. Again, the conclusion is that theoretical uncertainty is completely dominated by QCD 
effects at zero momentum transfer. 
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Finally, let us define an effective APV weak-mixing angle by the following relation: 

sin^ ^APV = 4v(0) sin^ HK)- (27) 

For = 91.1875 GeV, = 150 GeV and a^{Ml) = 0.119 we obtain sin^^APV = 0.231601 
(0.232123), corresponding to perturbative (non-perturbative) treatment. 

Appendix: Taylor expansions Here we list all expansions that are needed in order to 
reproduce the OMS results. Note that we need at most terms of O (s): 
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where M = M^,Mj^ and we remind a short hand notation for mass ratios: rij — — |- . 
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